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Abstract 

Coupled, multi-field models of inflation can provide several attractive features unavailable in the 
case of a single inflaton field. These models have a rich dynamical structure resulting from the 
interaction of the fields and their associated fluctuations. We present a formalism to study the 
nonequilibrium dynamics of coupled scalar fields. This formalism solves the problem of renormal- 
izing interacting models in a transparent way using dimensional regularization. The evolution is 
generated by a renormalized effective Lagrangian which incorporates the dynamics of the mean 
fields and their associated fluctuations at one-loop order. We apply our method to two problems 
of physical interest: (i) a simple two-field model which exemplifies applications to reheating in 
inflation, and (ii) a supersymmetric hybrid inflation model. This second case is interesting because 
inflation terminates via a smooth phase transition which gives rise to a spinodal instability in one 
of the fields. We study the evolution of the zero mode of the fields and the energy density transfer 
to the fluctuations from the mean fields. We conclude that back reaction effects can be significant 
over a wide parameter range. In particular for the supersymmetric hybrid model we find that 
particle production can be suppressed due to these effects. 

PACS numbers: 98.80.Cq, 11.15-z 



1 



I. INTRODUCTION 



In recent years, the study of nonequilibrium dynamics in quantum field theory has re- 
ceived much attention in various areas of physics, and particularly in cosmology. The work 
has been driven largely by inflation the most successful known mechanism for explaining 
the large-scale homogeneity and isotropy of the universe and the small-scale inhomogeneity 
and anisotropy of the universe 0. With observations for the first time able to directly test 
the more detailed predictions of specific inflationary models, the efforts in understanding 
inflation and its dynamics have redoubled. 

One area of particular interest is the dynamics of multi-field models of inflation in which 
the inflaton is coupled to another dynamical field during inflation. These models can lead 
to a variety of features unavailable in the case of a single field. Such multi-field scenarios 
include the well known hybrid inflation models 

On top of the dynamics during inflation, the subsequent phase of energy transfer between 
the inflaton and other degrees of freedom leading to the standard picture of Big Bang Cos- 
mology has been the subject of intense study. The inflaton may decay through perturbative 
processes 0, ^ as well as non-perturbative parametric amplification ||^, 0. The latter can 
lead to explosive particle production and very efficient reheating of the universe. 

Hybrid inflation and reheating models share an important common thread. They both 
involve the coupling of two or more dynamical, interacting scalar fields (or higher spin fields 
P]). An important aspect of such systems is the possibility of mixing between the fields. 
In Ref. H for example the classical inflaton decay is investigated for a two field model by 



solving the non-linear equations of motions on a grid. In Ref. |]T0|], the authors treat the 
problem of coupled quantum scalar and fermion fields at the tree level. Due to the small 
couplings involved in inflationary cosmology, such a tree level analysis is useful in a variety 
of physical situations. 

However, hybrid models as well as the dynamics of reheating typically include processes 
such as spinodal decomposition |]TT], O] and parametric amplification which require one to 



go beyond the tree level by including quantum effects either in a perturbative expansion or 
by means of non-perturbative mean field techniques such as the Hartree approximation or 
a large- iV expansion [|, [l^ . 



Going beyond tree level brings in the issue of renormalization. The problem of renormal- 
ization of time evolution equations in single field models was understood several years ago. 
In one of the first papers in this field. Cooper and Mottola showed in 1987 (Ref. [0, that 
it is possible to find a renormalization procedure which leads to counter terms independent 
of time and initial conditions of the mean field. They used a WKB expansion in order to 
extract the divergences of the theory. In a later paper Cooper et al. also discussed a closely 
related adiabatic method in order to renormalize the 0^-theory in the large N approxima- 



tion. Also Boyanovsky and de Vega, Ref. [11|, used a WKB method in order to renormalize 



time dependent equations in one-loop order, later on Boyanovsky et al. |T2[ investigated a 
(j)^ model in the large-N approximation and the Hartree approximation, too. In 1996 Baacke 
et al., Ref. |16[, proposed a slightly different method in order to extract the divergences of 
the theory, which enabled them to use dimensional regularization. In contrast to the WKB 
ansatz this method can be extended for coupled system, which was demonstrated in Ref. 



T7| . This procedure will be used also in this paper. We work in the context of a closed time 
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path formalism appropriate to following the time-dependent evolution of the system. 
In this formalism, the m-vacuum plays a predominant role, as quantities are tracked by 
their in-in expectation values (in contrast to the m-out formalism of scattering theory). We 
construct the m-vacuum by diagonalizing the mass matrix of the system at the initial time 
t = 0. However, because of the time- dependent mixing, a system initially diagonalized in 
this way will generally not be diagonalized at later times. 



One approach to this problem, taken in Ref. |]T0[, is to diagonalize the mass matrix at 



each moment in time through the use of a time-dependent rotation matrix. The cost of doing 
so is the appearance of time derivatives of the rotation matrix into the kinetic operators of 
the theory. While such a scheme is in principle workable beyond the tree level, the modified 
kinetic operators introduce complications into the extraction of the fluctuation corrections 
as well as the divergences that are to be removed via renormalization. 

We take an alternative approach where the mass matrix is allowed to be non-diagonal for 
all times t > and account for the mixing by expanding each of the fields in terms of all of 
the m-state creation and annihilation operators. The cost of doing so in an A^-field system 
is the need to track N"^ complex mode functions representing the fields instead of the usual 
A^. However, this allows standard techniques to be used to properly renormalize the system. 
For the two-field systems common in inflationary models, this effective doubling of the field 
content adds a relatively minor cost. 

For simplicity and clarity, we will work in Minkowski space time and in a one-loop ap- 
proximation. Extensions both to Friedmann- Robert son- Walker spacetimes and to simple 
non-perturbative schemes such as the Hartree approximation, while more complicated than 
the present analysis, present no fundamental difficulties. We note that Minkowski space 
time is a good approximation in the latter stages of certain hybrid inflation models, and 
it will also allow comparison with much of the original reheating literature |^, ^ which 
often neglects the effects of expansion, allowing us to directly determine the role played by 
the mixing of the flelds in the dynamics. 

The outline of the paper is as follows. We begin by considering the Lagrangian for 
coupled scalar flelds and set up our formalism for the quantization of the system. This is 
followed by an outline of the renormalization procedure. We then provide a summary of the 
results for the two-fleld case. We demonstrate the formalism with two examples: a simple 
reheating model and a hybrid inflation model motivated by supersymmetry. 

In the reheating model we investigate two relevant regimes discussed in detail in the 
literature p, 0], viz., the narrow resonance regime and the broad resonance regime. These 
different regimes occur depending on the choice of initial conditions. Usually in these models 
the mixing effects of the flelds were neglected by choosing a vanishing initial value for one 
of the mean flelds: We are now able to treat the full system and to investigate these mixing 
effects. For this purpose we concentrate on studying the behavior of the fluctuation integrals 
for the different flelds and the time-dependent mixing angle. Depending on the regime, as 
the mean flelds evolve, the effects of the mixing can be quite different. In the narrow 
resonance regime the mixing angle is very small and plays a sub-dominant role, whereas in 
the broad resonance regime the mixing effects are very important. Therefore, we emphasize 
that neglecting the mixing could lead to incomplete results. 

Supersymmetric hybrid models are a special realization of general hybrid inflationary 
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models (see e.g. Refs. [0, 0). Based on a softly broken supersymmetry potential, the 



special feature of these models is the occurrence of only one coupling constant, whereas in 
nonsupersymmetric hybrid models there are at least two different couplings. Thus, in the 
super symmetric case there is only one natural frequency of oscillation for the mean fields 
as long as fluctuations are neglected. This leads to efficient particle production during the 
preheating stage in the early universe. However, we show below that, by taking into account 
the fluctuations and investigating the full mixed system, this feature of supersymmetric 
hybrid models can be lost in some regimes. This is because the effective mass corrections 
for the two fields are different in these regimes, which leads to a chaotic trajectory for the 
renormalized field equations of motion in a phase space which mimics the situation of a 
nonsupersymmetric hybrid model. It appears, then, that supersymmetric hybrid models 
can lose some of their attractiveness compared to general hybrid models. 



II. N FIELDS 

We work with the following Lagrangian for real scalar fields $j with i = 1 . . . N: 

N ^ 

1=1 ^ 



where the potential is 



^11 1 

V{^i) = ^i^i + T^^ij'^i'^j + ■:Ti9ijk'^^'^j'^k + - A,jfei$,$j$fc$, . (2.2) 

i,j,k,l=l ^ ^■ 

Note that rriij, Qij^, and Xijki are symmetric in each index, but are generally non-diagonal 
resulting in the mixing of the different fields. In what follows, subscripts and superscripts 
i . . .n run over the values 1 . . . and we use a convention in which summation is assumed 
over repeated lowered indices, but not raised indices. 

We will expand each field about their expectation values (taken to be space translation 
invariant): 

<!>i{x) = Ut) + , Mt) = (H^)) ■ (2.3) 

Expanding the equations of motion and keeping terms to quadratic order in the fluctua- 
tions yields a one loop approximation. The equations of motion for the zero modes (pi are 
determined via the tadpole condition. We have 

4>i + Ai + niijcpj + ^gijk4>j(pk + ^Kjki4>j4>k4>i 
2 o 

+ ^^ijk^^j^k) + ^\jki4>j{S<l)kS4>i) = . (2.4) 
To this order, the fluctuations obey the equation 

6(f)^ - V^6(j), + M^j6(j)j = , (2.5) 

with the mass matrix ^ 

Mij = rriij + QijkCpk + -^>^ijki4>k4>i ■ (2.6) 
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As indicated in the introduction, the comphcation that arises is not the fact that the mass 
matrix (|2.6|) contains mixing between the various fields, rather that the mixing changes with 
time as the (pi evolve according to ( p.4| ). This means that if we diagonalize the mass matrix 
at one time, it will not generally be diagonal at any other time. 

Nonetheless, it is most convenient to quantize in terms of a diagonal system at the initial 
time t = 0. We define the matrix 

= O^kMuOl , (2.7) 

and the corresponding fiuctuation fields 

Xi = Oi,6(^, , (2.8) 
where Oij is an orthogonal rotation matrix. T>ij is diagonal at the initial time: 

Vi, = D'6ij , (2.9) 
without summation over the raised index i. The Xi obey the equations of motion 

Xi - V^Xi + VijXj = . (2. 10) 

We quantize the system by defining a set of creation and annihilation operators ajj(fc) 
and aci{k) where a = 1 . . . N corresponds to the m-state quanta of frequency 



UOr. 



VWTd^ . (2.11) 



As the mixing changes in time, each of the fields Xi is expanded in terms of all of the m-state 
operators. We have 

= Ey [a^ik)Unk,t)e^'-' + aUk)Urik,t)e-^'-^ . (2.12) 

The initial conditions for the A^^ complex mode functions are 

[/f (fc, 0) = 5f , f/f (fc, 0) = -zc^«f/f (fc, 0) . (2.13) 
It is convenient to define the fiuctuation integrals 

(X,X,) = E / 77rT-,V-Ur{k.t)Uf{k,t) , (2.14) 
a=i-' (27r)-^2u;ao 

from which it is straightforward to determine the contributions appearing in the zero mode 
equations (p.4|): 

{5<PM,) = OlOji{X,X,) . (2.15) 
It will also prove convenient to introduce the rotated couplings 

Gijk = gumOfjO'^^ , (2.16) 

AjjfcZ = KjmnO'^^O'^i , (2-17) 
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which allows us to write the zero mode equations as 

+ lGijk{XjXk) + ^Aijki<Pj{XkXi) = , (2.18) 

while the mode functions obey the equations 

Ut{k,t) + {k^ + V,,) Uf{k,t) = . (2.19) 

In addition to the equations of motion, it is useful to have an expression for the energy 
density of the system. This is particularly true when one completes numerical simulations of 
the system, since energy conservation is a powerful check of the accuracy of the simulations. 
After once again decomposing the fields into their expectation values and fluctuations, the 
energy density to one loop order is 

^ = + ^i4>i + ^mj4>i(pj + ^gijk(pi4>j4>k + ■^Kjki(f>i(l>Ak(f>i 
+ 1 (xf) + ^ ((vx,)') + \v,,{x^x,) , (2.20) 

where we've defined the integrals 

= f 4^s7^\Unk,t)\\ (2.21) 



a=l 
N 



a=l 



(277)3 2CU,0 

d^k e 
(27r)3 2cu,o 



T.^,^\Unk,t)\'- (2.22) 



A. Divergence structure and renormalization 



The mode integrals in the equation of motion defined by Eq. ( p.l4|) and in the energy 



density defined by Eqs. (|2.21|) , (|2.22|) are divergent and have t be regulated, allowing for 



a renormalization of the theory. We require a method of extracting the divergent terms 
appearing in the mode integrals, a nontrivial task, since the mode equations vary in time 
and they are coupled. Our aim is now, to find counter terms, which are independent of the 
initial value of the mean fields in order to formulate a finite theory. The correct choice of the 
initial condition for the fiuctuations guarantees that the theory is renormalizable. One way 
to extract the divergences of the mode integrals is due to a WKB method which allows for 
a high momentum expansion of the mode functions. However, when the fields are coupled, 
as in the present case, the usual formulation of the WKB expansion runs into difficulties 
which are yet to be resolved. 

An alternative method has been developed [|16], |l^, ^ which relies on a formal pertur- 



bative expansion in the effective masses and time derivatives of the masses of the fields. As 
such, it results in a series expansion of the mode functions in powers oim/uj and mjuP' ^ etc. 
The first few terms in the series include the divergent parts of the integrals that are to be 
removed via renormalization. 
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We begin by introducing the following ansatz for the mode functions: 

Ija ^ ^-iu^^ot ^^a ^ _ (2.23) 

The first term on the right hand side anticipates a quadratic divergence in the quantities 
{Xf). We define the following potential 

V^{t)=V,,{t)-D^6,,. (2.24) 

The equations of motion for the mode functions Eqs. (|2.19|) can be written in a suggestive 
form with the help of Eqs. (|]23H]24D 

/; - z2u;,o/; = - E (^r + fn ■ (2.25) 

1=1,2 

The terms on the right hand side of this expression are treated as perturbations to write the 
f's order by order in V, with the initial conditions /j(0) = /j(0) = 0. To first order in V, 
we have the equations of motion: 

ff^ - 2iu^off^ = -V^"„ . (2.26) 
The corresponding integral solutions for the real part of the /'s are: 

2gfJ/;(i) = + f dt'^^ cos(2.;„ot') , (2.27) 



while the imaginary part is of order and does not contribute to the divergences 
Using these results, we find quadratic and logarithmic divergences: 

(*-'f.)div = /(0(^%-4K5). (2.28) 

which must be removed via some renormalization procedure while also providing finite cor- 
rections to the parameters of the theory. 

To make the renormalization scheme explicit, we adopt dimensional regularization. We 
define the following divergent integrals 

"^'^ ^ - (2.29) 



(27r)3 2v/FT7? 167r 

/d?k 1 
7 — ^ ^ = I~M , (2.30) 
(27r)3 4(A;2 + ^2)3/2 31W, I ) 

where /i is an arbitrary renormalization point and /_3 carries the infinite contributions. In 
dimensional regularization I-sifi) is given by 



I-M = T^A- + ^^^'l \ ■ (2.31) 



The infinite part of (XiXj) is found to be simply 

(^i^i) infinite = -'Dijl-^^ix) . (2.32) 
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This leads to mass and coupling constant counter terms of the following form 



8Ai = 


2 


o Vr^/ .-7 J"- * 


6mij = 


i/. 

2 


-3(/^) [dikmgkmj + \jfci"^fc/] 


5gijk = 


h. 

2 


'?,{l^)9ilmMm,jk , 


5\jkl = 


2 


-si^f^^ ^ijmn^mnkl ■ 



(2.33) 
(2.34) 
(2.35) 
(2.36) 

It is important to notice, that these counterterms are independent of the initial conditions 
of the mean fields 0j. 

In addition to these counterterms, there are finite corrections of the parameters coming 
from the finite parts of the integrals ( |2.28D : 



(-^j-^jr)(;iiv finite 



1 



1 ^ , 

Vii In — . 



From this, we extract the following finite contributions to the couplings and mass: 



Arriij 



^9ijk 
^Kjkl 



1 



327r2 
1 

32^ 
1 



GijjD^ + gijkrrijk In — 



^ijkkD + {gikigkij + Kjkifnki) In — - 



87r2 
3 

32^2 



\ 1 El 



■ ^ijmn^mnkl 



In 



D" 



(2.37) 

(2.38) 
(2.39) 
(2.40) 
(2.41) 



These finite corrections are also contributing to the energy density. In addition we find 
a finite part due to the cosmological constant renormalization. 

The full, finite equations of motion become 

'^i + Ai + AAi + {niij + Amjj) + ^ {gijk + Agijk) (j)j(l)k + ^ {>^ijki + ^Kjki) (f^jMi 



+-(G..,+A.,,,0,)|:/^^ 



£0.^ CO. I ^a. COL I ret ^CK* i rciT/a 
Jfc + -^fe/i + (^l Ik + TT^^i 



. (2.42) 



B. Two Fields 

The two-field case is often encountered, and the physical applications we present in the 
next section are both in this category. It is therefore worthwhile to pause to look at a few 
details of such systems. 

We begin with a system of two real scalar fields $ and X 

^ = \{d,^? + \{d,Xf - V^($, X) , (2.43) 
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with the potential 



1 



4!' 



:^2^2 ^ -^2^2 ^ - ^4 ^ ^^4 ^ ^$2^2 



(2.44) 



This Lagrangian has the same form as that studied in the preceding section with the iden- 
tifications 



$1 


= 


$2 = 




mil 


= ml 


, m22 


= m 


■^1111 


= A, 


A2222 


= K , 




= 0, 


9ijk 


. 



^1122 = 9~ 1 A1112 — A1222 



(2.45) 



The remaining components of Xijki are determined by the fact that it is symmetric in each 
of its indices. 



The mass matrix Ai is 



M 



m2+A0V2 + ^7V/2 



9' 



(2.46) 



For two fields, the orthogonal rotation matrix can be written in terms of a single mixing 
angle 6. The matrix has the form 



O 



cos 6 sin 6 
— sin 6 cos 6 



(2.47) 



where the mixing angle is determined by the t = mass matrix Ai, Eq. ( p.46|) , through the 
relation 



tan^^ 



2-Mi2(0) 



-M22(0) - A^ii(O) + V(-M22(0) - -Mii(O))^ + 4A^f2(0) 



(2.48) 



The eigenvalues of O are the diagonal elements of the matrix V, Eq. ( |2.7| ), at the initial 
time: 



\ 

' 



with the values 



MniO) + A^22(0) + V(^22(0) - ^ii(O))' + 



A^ii(O) + A^22(0) - ^/{M22iO)-MniO)f + ^MhiO) 



(2.49) 

(2.50) 

(2.51) 



For general times, the mass matrix for the fields Xi and X2, writing cg = cos 6 and 
sg = sin 9, is 



clMii + 2s0CeMi2 + SIM22 seC0M22 - seCeMu + (c^ - sl)Mi2 V 
seCeM22- seCBMii + {cl- sl)Mi2 -2clM22 ~ ^seCeMu + sjMu }' ^ 
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The zero mode equations, before renormalization, read 

+ + + ^^^2 ^ ^ g^^.(^) ^^^^^.^ ^ Q ^ ^2.53) 

2 

X + + ^X' + + E = , (2.54) 

where 

\-^sgCe(l) + ^sgce(l) + ^{4- sj)x ^s1(l) + ^cj(f)- g^secex /' 

( Ph+^i4x + 9^secg(f) pecgx- ^■sgcgx + ^{c1 - Sl)(f)\ 
K^sgcex- ^sgcex + ^{cj- sl)(l) plx + ^slx- g'^secgcj) . ) 

The total energy density of the system including the fluctuations can be expressed as 



(2.55) 



C ■'-12, -'-•2,-'- 2/2,-'- 2 2,-^/4,'* 4, 5'^, 2 2 

+ \ (Xf) + ^ ((VX,)'\ + lv,,{X,X,) . (2.56) 



Now we have to formulate finite equations of motion and a finite energy density. We 
adopt the renormalization procedure of section [11 A| for the field case. By using the 
identifications ( ^.45| ) it is to derive the appropriate counterterms for the two field case from 
Eqs. ( p.33| )-( p736D . We find in particular 

= ^(^"^0 + 9^rnl)I^3{fx) , Sml = ^(c/^mj + Km\)I_^{^i) , (2.57) 

Of course we get also similar to Eqs. (|2.38|) - (|2.41| ) in the N field case finite corrections 
to the masses and couplings of the form 

AmJ = {d\XcI + ghl) + D^Xs] + g\l) + g^mlU + AmjL2} , (2.58) 

^rnl = {d\ksI + g^cl) + D^g^s] + kc]) + g^m^L2 + /tm^Li} , (2.59) 

^{A% + ,%}, (2.60) 

^ [9' (L, + L,) + \g'XL, + \g'^L,^ , (2.61) 

+ , (2.62) 



AA 
A/ 



where 



327r2 



Li = s^ln — + c^ln — , L2 = In — + In — . (2.63) 
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As a result of these finite corrections Eqs. ( p.58| - |2.6^ ), the total Lagrangian Eq. ( |2.43|) is 
also modified. This is exactly the renormalized Lagrangian which we needed. We also find 
an additional finite contribution to the classical Lagrangian given by 



AC 



D 



647r2 



4^1 _ D'') + In — [(A + + + + 2(mJ + mj) 



(2.64) 



and, the final zero mode equations for and x ^i-re given by 



dAC 



la 



+ frfr + 



(27r)3 2^,0 
S>k 1 ' 



Vn 



(27r)3 2cu,o 



Vr 



la 



(2.65) 



and. 



1 1 f\ f*. 

X + (m| + AniDx + -{k + Aft)x' + -{g^ + A(?2)0;^2 ^ __ 



(27r)3 2cu, 



la 



2u? 



0. 



(2.66) 



After writing down the finite zero mode equations of motion we also have to renormalize 
the energy density. Again by using the ansatz ( |2.23|) we can extract the divergent terms 
of the fluctuation integrals in Eq. ( |2.56| ). In addition to the quadratic and logarithmic 
divergences we find a quartic divergence. This leads to a counter term which acts as a 
cosmological constant and has the form 



5A = -(m^ + m^)J_3(m) . 

Altogether the divergent part of the energy density reads 
/_3(m) 



(2.67) 



div 



4 



{g'm'^ + \m'^W + {g^m^um'^x + (A + + ^g'Wx' 



(2.68) 



This expression leads of course to the same counter terms we found for the equations of 
motion, and therefore also to the same finite corrections to couplings and masses. Therefore 
it is straightforward to formulate a finite energy expression. 

Now, we are in a position to discuss the physical applications of our problem. This we 
shall do in the next section, but first, we introduce one more quantity that is convenient in 
discussing the degree to which the mixing plays a role in the dynamics. 
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C. Time-dependent mixing angle 



To better understand how the system evolves in time, it is useful to have a measure of 
how much the fields and x iiiix each moment, and how this mixing evolves with the 
system. To provide us with a measure of the mixing, we introduce a time-dependent mixing 
angle 6(t), which is defined in terms of the time-dependent mass matrix M.{t), Eq. ( |2.46| ): 



tane(t) 



2Mi2{t) 



M22{t) - A<n(t) + ^{M22{t)-Mii{t)f + AMUt) 



(2.69) 



III. PHYSICAL APPLICATIONS 

After setting up the technical framework, we are now in a position to investigate some 
relevant cosmological multi-field models for inflation. We begin our analysis with a simple 
two-field model often used for studying the phase of parametric amplification. (This phase 
occurs just after the completion of inflation in chaotic inflationary models |^.) This 
model provides a useful context to analyze the effects due to field mixing. Next we turn our 
attention to a supersymmetric hybrid inflationary model, which is of particular interest in 
cosmology. As discussed in the literature (see for example Ref. [0) particle production (and 



hence reheating) in these models is much more efficient compared to the nonsupersymmetric 
hybrid models. Until now the mixing effects in these models have not been treated fully 
including back reaction effects of the quantum fluctuations in the mean field approximation. 
This approximation does not take rescattering processes into account and therefore we can 
not address the problem of thermalization. 



A. Reheating 

The reheating phase in chaotic inflationary models is characterized by two different 
regimes, which depend on the chosen initial conditions: the first is the narrow resonance 
regime, while the second is the case of broad resonance. In order to investigate these regimes 
we examine two different parameter sets, where only the initial values for the zero modes 
are varied. We find significant differences in the behavior of the fluctuation integrals as well 
as the mixing angle in the two regimes. 

To make the analysis as simple as possible, we set A = k = as well as = 0. For the 
remaining parameters, we set = 1, which just acts as a unit of mass, and = 0.01. With 
these parameters, the case usually studied in the literature, x(0) = 0, does not introduce 
any mixing between the fields since the off-diagonal elements of the (p-x mass matrix are 
proportional to x- However, taking x(0) = may not always be the case. For instance, 
X field could take a large vacuum expectation value during inflation, provided x is treated 
as a field other than the inflaton. For the purpose of illustration we may consider a non- 
zero initial condition for x which is of order its effective mass g(j) at the beginning of the 
preheating stage and examine the consequences. The initial condition for in the narrow 
resonance regime is fixed by the condition g'^(j){0)/4m'^ = 0.1 (remember, g is fixed to be 
g = 0.001) and for the broad resonance regime by g'^(j){0)/4m'^ = 100. If we take 0(0) 
to be approximately the Planck scale as appropriate to the end of inflation, this would 
correspond to ~ and ~ lO^^GeV^ for the two respective cases. Note that 
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these parameters are chosen to depict the phenomena of interest during a time scale that 
can be accurately simulated. The results are, in any case, representative of two field mixing 
in the narrow and broad resonance regimes regardless of the precise parameter values in any 
particular model. 




FIG. 1: The log of the fluctuation integrals (Xf) 
(solid), {X1X2) (dashed), and {X^) (dotted) for 
the narrow resonance regime. 



FIG. 2: The log of the fluctuation integrals {XD 
(solid), {X1X2) (dashed), and (dotted) for 
the broad resonance regime. 



Fig. 1 shows the log of the three fluctuation integrals (Xf), {X1X2) and for the 

narrow resonance case. These are seen to be dominated by the exponential growth of 
while the other contributions grow more slowly. Therefore, the evolution is characterized 
by production of X2 particles. We now turn to the broad resonance regime, where things 
look quite different. Here, each of the fluctuation integrals grows rapidly as shown in Fig. 2. 
Significant mixing of the species occurs along with copious particle production. 
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FIG. 3: Time-dependent mixing angle Q{t) for 
the narrow resonance regime. 



FIG. 4: Time-dependent mixing angle Q{t) for 
the broad resonance regime. 



The behavior of the fluctuation integrals is consistent with the behavior of the time- 
dcpcndcnt mixing angle G. Here, the mixing plays a sub-dominant role in the narrow 
resonance regime (Fig.3) with the mixing angle remaining near zero. This means that 
X2 predominantly corresponds to the x held, such that the process is one of x particle 
production, which is as expected. Concentrating on the time-dependent mixing angle in the 
broad resonance regime. Fig. 4, significant mixing between the fields is observed. The rapid 
variation in the mixing angle indicates that mixing between the fields plays a very important 
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role in the evolution of preheating in the broad resonance regime. The large influence this 
has on the behavior of the system is clear from the evolution of the zero mode components 
in Fig. 5 and x{t) Fig- 6. 




FIG. 5: Zero mode evolution with fluctuation 
(solid line) and without fluctuation (dotted line) 
for (p{t) for the broad resonance regime. 



FIG. 6: Zero mode evolution with fluctuation 
(solid line) and without fluctuation (dotted line) 
for xit) for the broad resonance regime. 



B. Supersymmetric Hybrid Inflation 



We now consider a hybrid inflationary model where the flnite coupling of two flelds plays 
an interesting role in the termination of slow roll inflation ||^. The particular model we 
study is based on softly broken supersymmetry [ff9l with the potential 



V{a,N) = lmW + l.l{N' 



(3.1) 



The fleld a plays the role of an inflaton during inflation while the fleld is trapped in a 
false vacuum (A^) = 0. The inflaton rolls down the potential along the a direction to reach 
a critical value a = ac- Once a reaches its critical value, the effective squared mass for the 
N fleld becomes negative and consequently it rolls down from the false vacuum to its global 
minimum through the mechanism of spinodal instability [|ll|, |12|. Thus, inflation comes to 
an end and both the flelds begin oscillations around their respective minima given by 



0, 



N = V2i 



(3.2) 



This is the onset of the preheating stage, which has been discussed in the literature |T^, |2^, 
3[. The difference between the supersymmetric hybrid potential and non-supersymmetric 
hybrid potentials lies in different coupling constants. In Eq. ( p.l[ ), there is only single 
coupling parameter Kg, while in the non-supersymmetric version there can be at least two 
different coupling constants. The above potential, except for the a mass term, can be derived 
very easily from the superpotential for F-term spontaneously supersymmetry breaking: 



W 



fa {n^ - at) . 



(3.3) 



The appearance of a mass term for a in Eq. (|3.1|) is due to the presence of soft supersym- 
metry breaking. Its presence is essential for slow roll inflation to produce adequate density 
perturbations and also to provide a correct tilt in the power spectrum [|^. The height of 
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the potential during inflation is given by KgCxf . Similar potentials to Eq. ( |3.1j ) can also be 
derived from D-term supersymmetry breaking as discussed in Refs. [Q. In these models 
the critical value Oc and the height of the potential energy are related to the Fayet-Ilhpoulus 
term coming from an anomalous U(l) symmetry. As in any inflationary model, hybrid 
inflation is constrained by COBE This imposes the bound 

KsCXc ^ 1.27 X lO^^lr^l GeV , (3.4) 

where ri is one of the slow roll parameters which determines the slope of the power spectral 
index [ESI. For our purpose we flx it to be |?7| ~ 0.01. 



In order to discuss the details of the physics we mention here the equivalence between 
Eq. ( p.l|) and Eq. (p.2|) . This helps us to establish direct relationship with our earlier analysis: 

A =0, K = 6Kg , g = 2ks , = nia- , m^ = —2K^a^. (3.5) 

Notice, that is negative. An interesting feature of the hybrid model is that irrespective 
of the values of the parameters Kg, cXc, and m„, as long as they satisfy the COBE constraints. 



the behavior of the mean flelds follow a common pattern once they begin to oscillate [19 
First of all, the mass term for the a fleld, rricr, becomes less dominant compared to the 
effective frequency for the two flelds, which is given by the effective mass for the two flelds 
during oscillations 

rricsa = rneSN = ^K^ac > . (3.6) 

Hence, there is a single natural frequency of oscillation, thanks to supersymmetry. Since 
the masses of the flelds are the same at the global minima, there exists a particular solution 
of the equations of motion for the a and flelds. Their trajectory follows a straight line 
towards their global minima: 

N = ±V2{ac-a) . (3.7) 

The maximum amplitude attained by the a fleld is ~ ctc, while the other fleld attains 
a larger amplitude A^ = \/2ac- (We remind readers that a = corresponds to the point 
where the effective mass for A^ fleld changes its sign.) This is the point of instability which 
we need to discuss here. From Eq. (|3.1| ), we notice that prior to the oscillations of the flelds, 
and during the oscillations, the effective mass square for a is always positive. However, this 
is not the case for A^, and its mass square can be positive as well as negative even during 
the oscillations of the flelds, provided the amplitudes are large enough. If the amplitudes 
for (J and A^ are such that they satisfy Eq. ( |3.7|) , then the effective mass square for the A" 
fleld is in fact always negative for a < Oc- If the amplitudes are large enough such that after 
the second order phase transition the initial amplitude for a ~ cTc, it is then quite possible 
that near the critical point the effective mass for the fleld A^ vanishes completely. As far 
as the motion of the mean fleld without including fluctuations is concerned this does not 
provide any new insight. However, if the flelds are quantized then the perturbations in the 
fleld, especially for bN ^ grow exponentially because uj]^ in Eq. ( |2.26| ) becomes negative for 
sufficiently small momentum k. This shows that the vacuum is unstable near the critical 
point Uc- 
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Another intuitive way to appreciate this point is to consider the adiabatic condition for 
the vacuum. The adiabatic evolution for the zero mode evolution for N field is given by 
|ci;Ar| -C uj%. This condition is maximally violated at the point where the effective mass square 
for becomes zero, and, violation in adiabatic evolution of the zero mode for suggests 
that many fluctuations of 6N are produced during the finite period when the adiabaticity is 
broken [^. This explanation is quite naive because the overall production of particles and 



fluctuations depends also upon the global behavior of the zero mode fields. The effect of 
corrections due to fluctuations might affect the production of particles and this is the point 
we are going to emphasize in our numerical simulations. 

In some sense the hybrid model is quite different from chaotic inflationary models. In 
chaotic models, the inflaton field rolls down with an amplitude ~ l/(mt), where m is the 
mass of the oscillating field. However, in the hybrid model the amplitude of the oscillations 
die down very slowly, allowing many oscillations of the a and A^ fields in one Hubble time. 
Thus, one could expect large amplitude oscillations of the fields for a long time. This 
crucially depends on the parameter cTc. If o"c -C Mp, then we notice that effective masses for 
a and A^ fields during oscillations are much larger than the Hubble parameter. The Hubble 
parameter is given hj H ^ na'^/Mp during inflation, so, 

^eSa ^ rrieSN ^ ^ I (3 8) 

H H ac 

provided the scale of is quite small compared to the Planck mass, we can effectively 
neglect the expansion of the Universe. 

In the supersymmetric hybrid model there are two regimes of interest. Just after the mass 
square of the A^ field becomes negative, the fields begin to oscillate with an amplitude which 
decreases as oc l/t^. When the field amplitude drops below \N{t)/y/2ac - 1| < 1/3, the 
amplitude of the oscillations decreases as oc 1/t. In this regime, when the expansion of the 
Universe is neglected, the amplitude of the oscillations remains constant and the oscillations 
are harmonic: 

N(t) 1 , , , , 

-1 + 3 cos(meff.t) . (3.9) 

The corresponding evolution equation for the a field can be found from Eq. 

In this paper we are neglecting the expansion of the Universe. We will to concentrate upon 
two regimes: One with large amplitude oscillations which leads to the following parameters 

A = K = 24, g = 4ml = 0, m% = -16 x 10~^^ 
a(0) = 1.4 X 10-^ A^(0) = 1 X 10"^^ , (3.10) 

and, the other with small amplitude oscillations with the parameters 

A = 0, K = 24 X 10-^ ^ = 4 X 10-3 ml = 0, = -4 x 10"^^ 
a(0) = 0.24 X 10-^ N{0) = 0.66 x 10-^ (3.11) 

The coupling constants are dimensionless while the other dimensionful parameters denoted 
in Planck units. We find below a marked difference in the zero mode behavior of the fields 
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a and in these two cases depending on whether the fluctuations are taken into account 
or neglected. 

In parameter set p.lOp , we study the features of the fields with a large amplitude. This 
can happen when the a and fields begin their oscillations just after the end of inflation. 
As mentioned earlier, after the end of inflation the maximum amplitude attained by the 
mean fields can be quite large a = and A^ = v^o"c. This is precisely the initial condition 
we have chosen for the mean fields for our numerics, as shown in Fig. 7. The values for 
and (Tc can be evaluated from Eq. (PTH), which yields 



CTc. 



V2 X 10" 



(3.12) 
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FIG. 7: Evolution of cr (solid line) and N (dotted line) without fluctuations for parameter 
set (|3lO|). 
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FIG. 8: Trajectory for the fields a and N for 



parameter set ( 3.10 ) without fluctuations. 



FIG. 9: Trajectory for the fields a and for 
parameter set (|3.10|) with fluctuations. 



We notice that the evolution for a and A^ fields without taking into account the fluctua- 
tions are anharmonic, see Fig. 7, and, their trajectories in the a — N plane is a straight line, 
as shown in Fig. 8. However, switching on the fluctuations leads to a completely chaotic tra- 
jectory as shown in Fig. 9. The departure from the straight line trajectory is quite significant 
and it tells us that the renormalized zero mode equations have different contributions to the 
parameter 6^^ and to the effective mass of the A^ field. This mismatch in the frequencies of 
the zero mode equations for a and A^ leads to an irregular trajectory. 

The other way to interpret this behavior is to think in terms of different effective mass 
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corrections to cr and fields, such that the effective frequencies of the oscillations for a and 
do not match each other at the bottom of the potential. This is certainly a nontrivial 
result. Nonetheless, the result is quite expected from the fact that the amplitude of the 
oscillations are quite large and the effective mass for the N field is zero at each and every 
oscillation when a = ar- 
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FIG. 10: Zero mode evolution for a{t) with fluc- 
tuations for parameter set ( ^.10 ). 



FIG. 11: Zero mode evolution for N{t) with fluc- 
tuations for parameter set ( 3.10| ). 



As we mentioned earlier, the frequencies of the oscillations of the zero modes are different, 
as can be noticed in Figs. 10 and 11. The zero mode of a influenced by the fluctuations 
oscillates around its minimum a = with a more rapid frequency than when fluctuations are 
neglected. This suggests that the effective mass correction to the zero mode for a is coming 
solely from the flnite coupling contribution from the A^ fleld. (Note, that we have already 
set the bare mass for mg- = 0.) The oscillations maintain the regularity with increasing and 
decreasing amplitude. However, the story is not the same for the zero mode behavior for 
the A^ fleld. The amplitude of A^ increases gradually and the frequency of the oscillations 
varies. We mention here that the effective mass for the A^ fleld can vanish at a critical point. 
As a result, the adiabatic condition for the A^ fleld is violated at those instants and this is 
the reason why the amplitude of the A^ fleld is enhanced rather than suppressed. 

The evolution of the energy density is shown in Fig. 12. At flrst instance it seems quite 
odd that the energy density of the fluctuations does not increase further. One would naively 
expect a larger contribution of the energy density of 6a and 6N. This is not the case here. 
The energy density for the mean flelds and the fluctuations are equally shared. The reason 
is the correction due to the fluctuations. These corrections modify the effective mass of the 
A^ fleld and induce corrections to the coupling constants, namely k and g. The coupling 
constants are modifled in such a way that the trajectory of zero mode flelds become irregular. 
Usually the production of fluctuations is not efficient in this case. This is quite similar to 
the situation of preheating in non-supersymmetric hybrid models [22|. Even though we 



started with a supersymmetric hybrid model where at the bottom of the potential there is a 
single effective frequency, the situation changes completely if the fluctuations are taken into 
account. Essentially the coupling constants get a large correction which does not preserve an 
effective single coupling constant for the evolution of the zero mode flelds. This is precisely 
the reason why the zero mode trajectory becomes irregular and also why the production of 
66 and 6N is so low. 
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FIG. 12: Energy density stored in 5a and 5N fluctuations (dashed line) and zero mode 
energy density (dotted line) for parameter set ( p. 101) . The total energy remains constant 
(solid line). 

As a next example of the supersymmetric hybrid model we choose parameter set ( p. Ill ), 
with a small coupling and small dc 

Ks = 2xl0"^ = 0.707 X 10"^ (3.13) 

In this example, the coupling between the fields is quite small; g = 4 x 10~^, and also the 
initial conditions for a and have been chosen such that the fields oscillate around their 
respective minima. The maximum amplitude for a(0) = crj3 and A^(0) = (2v/2/3)(t, is much 
below the critical point cxc. We remind the readers that the chosen initial conditions for the 
oscillations do not come naturally just after the end of infiation, therefore this example 
does not represent a real situation. In spite of this we study the particular situation in 
order to notice the contrast in the behavior of the zero modes and the energy densities in 
the fiuctuations. This particular set of initial conditions for o"(0), A^(0) offers an alternative 
example where spinodal instability in the N field does not take place. As a result the 
effective mass square for the A^ field never crosses zero and the adiabatic condition for the 
A^ field is not strongly broken. The oscillations of the mean fields a and A^ are harmonic in 
nature, as shown in Figs. 13 and 14 by the dotted lines. The amplitudes are constant with 
a frequency given by Eq. ( p.6|) . The oscillations of the mean fields is governed by Eqs. ( p.7|) 

and (|3.9| ). The trajectory in the a A^ plane is a straight line whose slope is governed by 

Eq. (13). 

The effect of the fluctuations is also quite expected in this case. The amplitudes of 
the zero mode for a and A^ flelds decreases after a while and, in contrast to the preceding 
example, the frequency of the oscillations do not change very dramatically; see the behavior 
of zero mode in solid lines in Figs. 13 and 14. The trajectories for the zero mode evolution 
remain a straight line in this shown in Fig. 15. This is quite reasonable for the 

parameters we have chosen, but an important observation is that the effect of fluctuations 
does not alter the straight line trajectory for the zero mode flelds. This suggests that for 
small amplitude oscillations the corrections to the coupling constants; Ak, and Ag are such 
that the zero mode equations still have a similar oscillating frequency. This can be seen in 
Figs. 13 and 14. The production of Sa and ^A^ is not very signiflcant because the energy 
density stored in Sa and 6N does not grow rapidly. Thus the energy transfer from the zero 
modes to the fluctuation modes is not favorable for such a small amplitude oscillations as 
can be seen in Fig. 16. 
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FIG. 13: Zero mode evolution for a{t) with fluc- 
tuations (solid lines) and the without fluctua- 



tions (dotted line) for parameter set (3.11) 




FIG. 14: Zero mode evolution for x(i) with fluc- 
tuations (solid hues) and without fluctuations 



(dotted hne) for parameter set ( 3.11 ) 
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FIG. 15: Zero mode trajectories with fluctua- 
tions for parameter set ( |3.11| ). 



FIG. 16: Energy density stored in 6a and 5N 
(dotted hne) along with the zero mode energy 
density (double dotted line) for parameter set 4. 
The total energy density remains constant (solid 
line). 



We conclude this section by mentioning that preheating in this supersymmetric hybrid 
model is quite interesting. Depending on the amplitude of the oscillations of the fields, the 
behavior of the zero mode can be quite different. As a new feature we noticed that if the 
amplitude of the oscillations is close to the critical value Uc, the effective mass square for the 

field becomes negative and as a result the fluctuations of the field grows exponentially. 
However, the effect of fluctuations alters the coupling constants in such a way that the 
trajectory of the zero modes become irregular. Even though the adiabatic conditions seem 
to be broken for the field near the critical value, the energy density transferred from the 
zero mode to the fluctuations is not sufficient. Our study reveals some interesting messages 
which we briefly mention here. We emphasize the point that the departure from the straight 
line trajectory of the zero mode is an essential feature of a supersymmetric hybrid model 
if the fluctuations are taken into account. Even though, we have not included the Hubble 
expansion, the results we have obtained are quite robust because supersymmetric hybrid 
inflationary models have a unique behavior of the fields which allows a smaller inflationary 
scale compared to the effective masses of the fields around their global minima. This suggests 
that during the oscillations, the expansion is felt much later, on a time scale determined 
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by the parameters. This behavior is not shared by models where inflation is governed by a 
single field as in chaotic infiationary models. This undermines the production of quanta from 
the vacuum fluctuations. In several ways this affects the post inflationary radiation era of 
the Universe. Supersymmetric, weakly interacting dark matter formation and generation of 
baryonic asymmetry in the Universe during preheating are the two most important frontiers 
which due to our results may warrant a careful revaluation. 

In order to substantiate our claim that a due consideration of fluctuations after the end 
of inflation is an important feature of any supersymmetric hybrid model, we have chosen 
an unphysical example which serves the purpose of making a vivid distinction. We stress 
here that the spinodal instability which is actually responsible for producing an irregular 
trajectory of the zero mode of the fields in a phase space is completely lacking if the ampli- 
tudes of the oscillations for o", are small compared to the critical value a^.. This acts as 
a comparative study and shows that after the end of infiation, in a supersymmetric hybrid 
infiationary model, due to the spinodal instability in a field, a proper renormalization of the 
masses and the coupling constant have to be taken into account. 

IV. CONCLUSION 

We have introduced a formalism to address the dynamics of N nonequilibrium, coupled, 
time varying scalar fields. We have shown that the one-loop corrections to the mean field 
evolution can be renormalized by dimensional regularization. For the sake of clarity and 
simplicity we restricted ourselves to Minkowski spacetime while deriving the renormalized 
equations of motion and the energy density of the system. We applied our formalism to a 
two field case where we study the behavior of the quantized mode functions and the effect 
of ffuctuations on the zero mode equations of motion for various parameters including small 
and large amplitude oscillations and large and weak coupling between two scalar fields. The 
varied couplings and amplitudes illustrate various facets of the intertwined dynamics of the 
two fields which lead to a deeper understanding of the production of self quanta and transfer 
of energy density between the fields in a cosmological context. 

As a special example we have chosen a two field infiationary model which is genuinely 
motivated by supersymmetry and thus preserves the effective masses of the fields to be the 
same in their local minima. The model, as a paradigm, predicts inflation which comes to an 
end via a smooth phase transition, and robustness of the model is conflrmed by a slightly 
tilted spectrum of scalar density fluctuations within the COBE limit. The model parameters 
can be adjusted to give an inflationary scale covering a wide range of energy scales from TeV 
to 10^^ GeV. The phase transition leads to a spinodal instabihty in one of the fields which 
leads to a coherent oscillations of the fields around their global minima. The instabihty 
occurs in one of the fields which demands careful study of the back reaction to an otherwise 
growing mean field in an intertwined coupled bosonic system. An account of infiuence of 
the fiuctuations gives rise to uneven contribution to the renormalized masses of the fields. 
This results in an irregular trajectory of the zero mode in a phase space, which breaks the 
coherent oscillations of the two fields. This prohibits an excessive production of particles 
from the vacuum fiuctuations. This requires a careful revaluation of the successes of the 
production of weakly interacting massive particles and baryogenesis via out of equilibrium 
decay in supersymmetric hybrid inflationary models. Our study implies that exciting higher 
spin particles from the vacuum fluctuations of the coherent oscillations of the fields in a 
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supersymmetric hybrid inflationary model demands careful reconsideration. 

Even though, we have neglected the effect of expansion in our calculation, our results are 

robust enough to claim that the fluctuations in a supersymmetric hybrid model do not grow 
if the back reaction of the fluctuations are taken into account in the mean field evolution. 
An extension of our formalism to an expanding Universe deserves separate attention. 
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